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ABSTRACT 

Recent analyses of the Wilkinson Microwave Anisotropy Probe (WMAP) data have suggested 
C") ■ that the low-order multipoles of the cosmic microwave background (CMB) anisotropy distri- 

bution show cosmologically interesting and unexpected morphologies and amplitudes. In this 
paper, we apply a power equalization (PE) filter to the high latitude WMAP data in order to 
reconstruct these low-^ multipoles free from the largest Galactic foreground modelling un- 
certainties in the Galactic plane. The characteristic spatial distributions of the modes of order 
• I — 2, 3,4, 5 have been determined as a function of frequency, sky c overage and two meth- 

ods for foreground correction: using the template-based corrections of iBennett et alJ (l2003bl) 
f^S . and by a simple linear projection scheme assuming the spectral dependence of the foreground 

components. Alth ough the derived multipole maps are statistically consiste nt with previou s 
estimates from Tesmark, de Oliveira-Costa & Hamilton (2003) (TOH) and Efstathiou (2003), 
our analyses suggest that the K and Ka frequency bands remain significantly contaminated by 
residual foreground emission for the WMAP Kp2 mask. However, the I = 3,4,5 multipole 
maps for the Q, V and W channels indicate that, after foreground cleaning, these multipoles 
are dominated by the CMB anisotropy component. We confirm the Teamark et al. (2003) re- 
sult that the octopole does indeed show structure in which its hot and cold spots are centred 
on a single plane in the sky, and show further that this is very stable with respect to the applied 
mask and foreground correction. The estimated quadrupole is much less stable showing non- 
r> | negligible dependence on the Galactic foreground correction. Including these uncertainties is 

likely to weaken the statistical significance of the claimed alignment between the quadrupole 
and octopole. Nevertheless, these anisotropy patterns are also present in the COBE-DMR data, 
and are unlikely to be associated with instrumental systematic artifacts. 

Key words: cosmic microwave background-cosmology: observations 



o 
o 
in 



Oh 



1 INTRODUCTION 

In the standard cosmological interpretation, the distribution and 
statistical properties of the cosmic microwave background (CMB) 
temperature anisotropies reflect the properties of the Universe ap- 
proximately 300 000 years after the big bang. The simplest in- 
flationary models of the origin of the primordial density fluctua- 
tions predict statistically isotropic and Gaussian fluctuations for the 
CMB. Thus, in the conventional approach, where the temperature 
anisotropy field, AT(9, (f>), is expanded as a sum over the spherical 
harmonic basis on the sphere 

AT(e,4>) = Y J ^] aemYemie '^ (1) 

1=1 m =-l 
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the decomposition coefficients, a[ m , are independent, Gaussian dis- 
tributed stochastic variables. Then, as a result of statistical isotropy, 
their variances Q = (\a {m \ 2 ^ depend only on t. The validity of these 
predictions can be constrained by direct observation of the CMB 
anisotropy. 

The sum over m for a given I, AT e (6, 0), is then a multi- 
pole map of order I, and typically represents the CMB fluctua- 
tions on angular scale ~ 1/f. The lowest multipoles of the CMB 
anisotropy, in particular the dipole (I = 1), quadrupole (I = 2) 
and octopole (I = 3), are especially important for studying the 
homogeneity and isotropy of the Universe afte r the last s catter- 
ing of the CMB photons iMaar tens. Ellis & StoegeJll996l) . The 
measurement of these low-order multipoles requires observations 
over a large fraction of the sky, as provided initially by the COBE- 
DMR experiment, and more recently by the Wilkinson Microwave 
Anisotropy Probe (WMAP) satellite. The dipole is generally inter- 
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preted as being dominated by a kinetic effect related to the peculiar 
motion of the Earth with respect to the CMB rest frame, so cos- 
mologically relevant studies usually consider multipoles only from 
the quadrupole upwards. The first measurement of the quadrupole 
by COBE-DMR suggested an intriguingly low value for its am- 
plitude, and l ed to vigorous discussion about the credibility of 
its estimation JSmoot et alJll992l : iBennett et ailll992t IStarldll993t 
iGould l 1993; Kogut et alJ ll996l) . In particular, doubts centred on 
limitations imposed by systematic effects, instrumental noise, and 
non-cosmological foreground contributions due to local astrophys- 
ical emission. Nevertheless, the low amplitude of the quadrupole 
was recently confirmed by high precision measurements from the 
WMAP satellite l Bennett et al. 2003a) after the application of a re- 
fined analysis of Galactic foreground contamination. 

Subsequently, iTeemark et alj 120031) (hereafter TOH) have 
pointed out that the quadrupole and octopole appear rather planar, 
with most of the hot and cold spots centred on a single plane in the 
sky. Moreover, the two planes appear rather aligned. Other stud- 
ies have also demonstrated the presence of intriguing correlations 
ICopi. H uterer & Starkmanl2004|:[Schwarz et al 12004 1) and anoma- 
lies jEriksen et alj l2004bt Larson & Wandeltl 120041 IVielva et"ai] 
120041) of the low-order WMAP multipoles. Some of these features 
could be explained by a specific non-standard topology for the 
Universe, although a broad class of finite universe models, includ- 
ing the Poincare Dodecahedron model proposed bvlLuminet et alj 
(2003), was recently ruled out bv lCornish et alj 120041) . However, 
many of these previous analyses have determined the multipole am- 
plitudes on the full sky anisot ropy as m easured by either the inter- 
nal linear c ombination (ILC. IBennett et al 12003^ ' or TOH map. 
Importantly, Eriksen et al. 1 2004a) have shown that residual Galac- 
tic foreground uncertainties for the ILC map in the Galactic plane 
can be at the level of a pproximately 50 yuK, whilst the analysis of 
ISlosar&Seliakl (2004) strongly suggests that much of the detected 
alignment is indeed related to that part of the data most contami- 
nated by (residual) Galactic foreground. 

The goal of this paper is to estimate the low-order WMAP mul- 
tipoles using a technique which allows the explicit determination 
of their reconstruction uncertainties, and subsequently to clarify 
whether the unexpected features of the quadrupole and octopole 
are significant. Some of these uncertainties are related to Galactic 
foreground emission. In what follows, we consider that the estima- 
tion of CMB anisotropy based on the WMAP maps is not credi- 
ble when including data from the Galactic plane, and we will omit 
this region using a set of increasingly aggressive masks. Even after 
this data excision, the remaining high-latitude Galactic foreground 
emission may compromise our multipole reconstructions. There- 
fore, we adopt two independent methods of foreground removal. In 
the first, we use three templates for the Galactic synchrotron, free- 
free and dust (including the so-called anomalous dust-correlated 
emission), scaled to the WMAP frequencies by the weights pro- 
vided by the WMAP team and subsequently subtracted from the in- 
dividual observation maps. This treatment is completely consis- 
tent with that of the WMAP team and has largely been adopted 
as the standard correction by the CMB community at-large. For 
the second method, we do not use templates, but simply take lin- 
ear combinations of the WMAP data in five frequency bands de- 
signed to project out the CMB component after assuming spectral 



Indeed, it should be noted that such analyses explicitly ignore the ad- 
monition of the WMAP team that the map is not suitable for cosmological 
analysis. 



index values for the individual foreground components. More de- 
tails about the methods and their limitations will be found in the 
appropriate sections of the paper. The multipoles determined from 
the foreground-corrected data will also be compared to those com- 
puted from the ILC and TOH cleaned maps, to which we also apply 
masking of the Galactic plane. Because these were derived using a 
different technique for removing the Galaxy, this provides yet fur- 
ther constraints on the foreground contributions. Note that it is not 
our intention to provide a definitive model of Galactic emission, but 
simply to investigate the stability of the inferred low-^ anisotropy 
patterns to different assumptions. If the multipoles are stable to 
foreground correction then we will infer that the anisotropy is likely 
to be cosmological in origin. 



2 THE PROBLEM 

One of the major goals of measurements of the CMB radiation is to 
allow for the construction of a map of the anisotropy. To facilitate 
this, the measurements are discretized into _/V p j x pixels, which can 
be written as a Af p j x -dimensional data vector d. This is related to the 
true sky signal x (itself the sum of the true CMB anisotropy field 
xcmb and the foreground contribution Xf g ) by 

y = Bx + n , (2) 

where the matrix B encodes the effect of beam and pixel smoothing, 
while the additional term n represents the contribution from detec- 
tor noise. The foreground contamination can initially be minimized 
by masking out those regions of the sky with strong foreground 
emission, e.g. as a result of the Galactic plane or point sources. The 
remainder of the map can then be corrected for the high-latitude 
Galactic (and extragalactic if necessary) emission. The vector d is 
then considered to denote a cut sky, foreground corrected map. 

In the spherical harmonic, Y{ m , space, the vector 2 of decom- 
position coefficients of the map on the cut sky, a, is related to the 
corresponding vector on the full sky, a, by 

a = K B a + n c , (3) 

where n r denotes the noise vector on the cut sky and B now con- 
sists of the appropriate harmonic coefficients of the beam window 
function. We will assume that the beam profile is azimuthally sym- 
metric. Then B is a diagonal matrix dependent only on £ and con- 
sisting of the weights wj! cam . To simplify notation, hereafter a will 
denote the vector multiplied by B. K is referred to as the coupling 
matrix and its (;, f)-th component is given by 

Ki(,ejnM?M) = I Y' Cm (n)Y [ , m ,(ri)dQ. il . (4) 

J cut sky 

Our aim is to determine the full sky representation of the 
anisotropy written in terms of the coefficients a using the above 
relation. Unfortunately, we do not know the noise contribution to 
this specific set of observations, only a statistical description of the 
ensemble properties, and moreover the sky cut causes the coupling 
matrix K to be singular. This means that there is no information 
on those anisotropies that project onto modes in the spherical har- 
monic basis where the only structure lies within the sky cut. Hence 
it is impossible to reconstruct all of the modes from the a. However, 
for low-order multipoles, small sky cuts and a high signal-to-noise 

2 Mapping from an index pair (t, m) into a single index i is given by 

i = e 2 + e + m + 1 
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ratio, a good approximation is to simply truncate the vectors and 
coupling matrix at some value of the index i(£, m) then reconstruct 
the multipole coefficients by inverting the non-singular truncated 
matri x K. This m ethod was applied to the first year WMAP data in 
lEfstafhioul iiool . and will hereafter be referred to as direct inver- 
sion. 

Direct inversion is a special case of a more general technique 
for estimation of the true sky signal by a linear transformation of 
the data vector a 



The subscripts £, and H denote the range of indices 



a = F a. 



(5) 



where matrix F acts as a filter and satisfies certain conditions im- 
posed on the solution a. In this paper, we will consider both lWiened 
( 1949) and power equal ization (PE) fil ters, and utilise the detailed 
formalism developed by Gorski 1 1997) for the application of these 
filters to CMB observations. 



2.1 Wiener filter 

The Wiener filter JWienej[l949T) is constructed from the require- 
ment that the mean square deviation between the filtered, a, and 
the true signal, a, e - ((a - a) r • (a - a)^ = tr((a - a) • (a - a) 7 ^, 
is minimized. Then we obtain 

F w = S ■ K r • (K • S • K r + n)~' , (6) 

(see, e.g jBunn et alll994l:lzaroubi et alll995l) where S = /a • a r ) 

and N = (n c • are covariance matrices of the signal and noise, 
respectively. One should notice that this filter practically anni- 
hilates those components of the data which are assessed to be 
noise dominated. Non-cosmological modes can be removed from 
the analysis by assumption that their power spectrum in the ma- 
trix (K ■ S • K T + N)~ is very large. In practice, the monopole and 
dipole are such terms, because for differencing experiments such 
as WMAP the monopole may correspond to an unphysical quantity 
related to the map-making process, whereas the dipole is attributed 
to the motion of the Earth in the Universe. For negligible noise the 
filter is essentially equivalent to the direct inversion method and 
does not depend on the assumed signal covariance matrix. How- 
ever, we would like to explicitly take into account the correlations 
induced by the cut between the lower order modes we are interested 
in reconstructing with the higher order modes present in the data. 
These correlations are non-negligible and play the role of noise in 
the Wiener filter. 

For further consideration, it will be useful to introduce a new 
orthonormal basis of functions on the cut sky. A method of con- 
structio n such of such orthonormal functions, i/j, was proposed by 
lG6rskilil99l . 

In this new basis, the vector of decomposition coefficients of 
the cut sky map, c, is related to vector a by 



'a + iu 



(7) 



where L is the matrix derived by the Choleski decomposition of the 
coupling matrix and is the vector of noise coefficients in the ip 
basis. Because L is a lower triangular matrix, c can be decomposed 
into a part dependent on both lower ax and higher a^ order modes 



c x = L T ££ ■ a £ + L T H£ ■ a w + n^x > 

and a part dependent only on higher order modes 



(8) 



(9) 



1, . . . , (4* + l) 2 and i = (4* + l) 2 + 1, 



+ 1) corre- 



sponding to those modes which are to be reconstructed with the 
filter (multipoles in the range from 1 to 4c) and the remainder 
of the modes used in the analysis (multipoles in the range from 
4c + 1 to 4iax)> respectively. It should be noticed that only the first 
part of the vector c contains information about lower order modes 
ax, so to estimate them we will use only this part of the vector. 
This is a significant difference as compared to the vector a where 
information about the lower order modes is also contained in the 
second part, a^. 

Taking Cx as a data vector, estimation of the lower order 
modes will be given by fix = F • C£. For the Wiener filter we have 



* w - &££ • ^££ • 
where C = L ££ 



(10) 



Lxx + L^ £ • Svh ■ L-«£ + and S ££ , 



S-hh are the appropriate covariance matrices of the CMB modes 
corresponding to the lower and higher multipoles, respectively. 
is the noise covariance matrix in the i/i-basis. 



2.2 Power equalization filter 

The PE filter is defined by the requirement that F is chosen such 
that 



a • a r ) = S££, i-e. F PE ■ C • F T PE = S ££ . 



(11) 



Using the Choleski decomposition of the relevant matrices, we can 
construct the PE filter as follows: 

si = l s1 ■ i£ v r = L -i , s £X = • r s , , ( i 2) 



Zl £1 



£€ ££ 



ZL ZZ 



and as a result we obtain an upper triangular matrix: 



J c-' 



(13) 



(14) 



where the matrices S ££ and C are the same as in the Wiener filter 
case. 

The condition II H means that on average over many ap- 
plications the PE filter renders filtered data whose statistically 
most likely power spectrum matches that of the underlying 
signal. Because the PE filter is an upper triangular matrix, the 
non-cosmological monopole and dipole terms can be explicitly 
eliminated from the analysis. Then, the restored modes (with ( > 4) 
are not contaminated by the monopole and dipole contributions. 



2.3 Errors 

The errors of estimation are determined on the basis of 
the diagonal terms of the covariance matrix, defined as 
M = ((ax - ar)(ax - fir) 7 /- Thus, we have 



M 



(F-L ££ -l).Sxx-(F-L X£ -l) r + 

+F ■ (IX, • S ww ■ Lwx) ■ F T + F • N„ ■ F r , 



(15) 



where I is the identity matrix. This is a sum of three components: 
the first is a result of correlations between the recovered low-order 
multipoles, the second is induced by correlations with higher order 
multipoles and the last is generated by noise. 
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2.4 General considerations 

In order to construct the filters we have to assume the form of 
both the signal and noise covariance matrices. We will make the 
usual assumption that both the CMB and noise components are 
Gaussian stochastic variables. The variances for individual modes 
of the CMB anisotropy, Q = are uniquely expressed as 

integrals over the power spectrum of the matter density fluctua- 
tions and depend only on C as a result of the statistical isotropy 
of the CMB temperature field. The rms noise level in pixel p, 
°"(p)noisc = o~ol V^otsCpX depends on quantities specific to a given 
instrument and scan strategy: the number of observations in pixel, 
NdbAp), and the rms noise per observation, cr . Dependence on 
the assumed power spectrum might seem to be a disadvantage of 
these filtering methods. However, we will see that in the case of the 
WMAP data, the multipole estimation does not depend significantly 
on the assumed power spectrum. 



3 DATA 

The WMAP satellite observes the sky in five frequency bands de- 
noted K, Ka, Q, V and W, centred on the frequencies of 22.8, 
33.0, 40.7, 60.8 and 93.5 GHz, respectivel y. The maps 3 are pix- 
elized in the HEALPix 4 scheme <G6rski. Hivon & WandelHl 19981) 
with a resolution parameter N s n c = 512, corresponding to 3 145 
728 pixels with a pixel size of ~ 7 arcmin. The maps are cor- 
rected for the Galactic foreground using three templat es 4 [a radio 
surve y at 408 MHz l Hasl arn et al Jl 198 lh . a Ho- map (Finkbeiner 
120031) . and a dust map based on the combined COB£-DIRBE and 
IRAS data l Finkbeine r. Davis & Schle gelll 19991)1 to t race th e syn- 
chrotron, free-free and dust emission. iBennett et alJ |2003b) have 
concluded that the combination of these templates, each scaled to 
the frequencies in question by a single factor derived for the high- 
latitude sky (i.e. making no allowances for spectral variations in the 
foreground behaviour), is sufficient to account for almost all of the 
foreground emission including the anomalous dust-correlated com- 
ponent (irrespective of its physical origin). Although the WMAP 
team did undertake an analysis of the data based on the maximum 
entropy method (MEM) which allowed for spectral index varia- 
tions of the foreground components on a pixel-by-pixel basis, the 
resultant CMB map was deemed inappropriate for use as a result of 
complex noise properties. For all cosmological data analysis, the 
WMAP team have corrected the data using the template method. 

We adopt the coefficients from Bennett et al. 1 2003b), which 
were derived by fitting the templates to the data for the Kp2 sky 
coverage, and then subtract the appropriately scaled data from the 
WMAP sky maps 5 . For the Q, V and W bands five coefficients were 
determined: an amplitude for each of the synchrotron and free-free 
components at 41 GHz (and scaled to the other frequencies by as- 
suming spectral indices of -2.7 and -2.15 respectively), and three 
for the dust, one at each frequency. The anomalous dust-correlated 

3 Available at http://lambda.gsfc.nasa.gov 

4 http://www.eso.org/science/healpix 

5 In fact, the LAMBDA web provides a set of foreground corrected maps 
generated with this recipe. However, a problem was discovered - see 
http://lambda.gsfc.nasa.gov/product/map/IMaps_cleaned.cfm - in that the 
templates used in this process were convolved with a too low /max (the 
maximum multipole component), resulting in ringing around strong point 
sources. Although we do not expect this to compromise our results on large 
angular scales, we have, nevertheless, generated our own template correc- 
tions. 



emission was thus effectively absorbed into the dust template scal- 
ing (which also accounts for the thermal contributi on). For the K 
and Ka channels, the coefficients from table 3 of Bennett et alJ 
(2003b) were used. The Q-, V- and W-band maps are then com- 
bined using inverse-noise- variance weights to create the corrected 
co-added map. 

We utilize several masks in the work which follows: the Kp2 
mask (eliminating ~ 15 per cent of the sky), the KpO mask (elim- 
inating ~ 23 per cent of the sky), an extended Galactic mask and 
mask 30. The extended Galactic mask (hereafter denoted 20+) is 
the same as that used for the COBE-DMR data analysis and is de- 
fined as the region with latitude \b\ > 20° plus custo m cut-outs 
in the vicinity of Orion and Ophiuchus (Bennett et al. 1996). The 
broadest mask 30 is defined as the region with latitude \b\ > 30°. 
The latter masks were additionally modified by excluding the point 
source mask included in the Kp2 and KpO masks, so that they fi- 
nally eliminate ~ 38 and ~ 50 per cent of the sky, respectively. 
It should be noted here that changing the Galactic cut applied to 
the data could, in principle, require modulation of the scaling co- 
efficients used for the templates, particularly if the foregrounds do 
demonstrate strong variations in spectral index with Galactic lati- 
tude. However, we proceed with the simple foreground model spec- 
ified by the WMAP team and used by them (and elsewhere in the 
literature) irrespective of the applied mask. 



4 CHOICE OF 4iax AND 4ec 

The choice of £ mix and £ IX is based on Monte Carlo (MC) sim- 
ulations of the WMAP data, consisting of both realizations of the 
CMB anisotropy and noise. The CMB anisotropy maps were gen- 
erated as random realizations of a Gaussian field with an angular 
power spectrum corresponding to the best-fitting ACDM cosmo- 
logical model with running spectral ind ex 6 iBennett et alj|20033 : 
ISpergel et all2003tlHinshaw et alj2003T) . convolved with the chan- 
nel specific beam window function wjp" 11 . The noise maps were 
produced from Gaussian random numbers with pixel-dependent 
variance corresponding to the properties of the corrected co-added 
map. The simulated maps were studied after masking by the Kp2 
mask. 

The upper limit for possible values of f max and t mc is im- 
posed by the feasibility of the Choleski decomposition of the cou- 
pling matrix. For a given cut, the coupling matrix becomes numer- 
ically singular for multipoles larger than fchoiesM at which point the 
Choleski decomposition fails. For the Kp2 mask this multipole is 
found at fchoieski ~ 70. The studies in this paper are focused mostly 
on estimation of the lowest multipoles (2 < t < 10), and partic- 
ularly on those modes which have been shown to have potentially 
anomalous properties (2 < t < 5). Thus it is sufficient to include 
only multipoles up to ^ max = 30. Increasing f m . ix to 35 or 40 did 
not appreciably change the estimated values of these multipoles. 
This is a direct consequence of the decreasing amplitude of higher 
multipoles, scaling roughly as £~ l (Q oc (~ 2 ). 

In order to establish a suitable value for 1000 MC simu- 
lations were generated and the PE filter computed with £ max = 30 
applied to the simulated maps masked by the Kp2 mask. The power 
spectrum of the difference between the true and restored multipole 



6 Available at http://lambda.gsfc.nasa.gov 
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coefficients over the range 2 < £ < 5 

yif /i OUtp 

c error _ f^'=- f 



"I 2 ) 



2f + 1 

is shown in Fig. Q as a function of £ rec . Notice that for 
Acc = fmax = 30, the errors are the largest. In this case, the PE fil- 
ter does not adequately take into account the aliasing of higher or- 
der modes into the range I < £ Iec . Significant errors are seen also 
for lower values of caused by correlations of the multipoles 
I = 2, . . . , 5 with modes of order l KC + 1, £ IK + 2 which for low 
values of £ mc have significant amplitudes. In what follows we adopt 
Acc = 10; although the flat power spectrum of the uncertainties over 
the range 10 < l mc < 25 allows a larger value for £ rx to be selected, 
later in the paper we will consider the broader Galactic cut 30 for 
which the PE filter fails when too large a value for f rec is used. 

It should be noted that other methods for the determination of 
the low multipoles also require the choice of some parameters to fa- 
cilitate the reconstruction method. In particular, the direct inversion 
method of Efstathiou (2003) adopts two parameters without justi- 
fication: the order of the multipole at which the coupling matrix is 
truncated and the FWHM of the beam for the initial smoothing of 
the data. 



5 COMPARISON OF THE WIENER AND PE FILTERS 

Fig. J3J provides a comparison of the WMAP power spectra ob- 
tained by the Wiener and PE filtering. In addition, we show the 
averages over 10 000 WMAP map simulations for the power spec- 
tra of the restored maps. In both cases the Kp2 mask and filters with 
f rec = 10, £ max = 30 were applied. 

The Wiener and PE filters differ mainly in the recovered power 
of the multipoles: there is a systematic suppression of power with 
increasing t for maps reconstructed by the Wiener filtering. Con- 
versely, the definition of the PE filter ensures an unbiased estimate 
for the multipole amplitudes. However, the PE filter does not give 
an optimal estimation in the minimum variance sense, i.e. an es- 
timate which has the least variance between the true and restored 
data, which is the case for the Wiener filter case. Nevertheless, we 
decided to use the PE filter in the remainder of the paper because 
for the lowest multipoles and Kp2 mask the results do not differ sig- 
nificantly from those derived by the optimal Wiener filter, and the 
unbiased nature of the estimator is considered to be more impor- 
tant. This is particularly relevant with increasingly aggressive sky 
cuts where the Wiener filter demonstrates even more suppression 
of power. 



6 COMPARISON OF THE DIRECT INVERSION AND 
THE PE FILTERING METHODS 

10 000 simulations of the WMAP data were again used to compare 
the accuracy of reconstruction of the low-order multipoles by the 
direct inversion and the PE filtering methods. The reconstructions 
were performed for the Kp2 mask, and for a PE filter constructed 
with ( teo = 10 and f max = 30. For the direct inversion method the 
coupling matrix was evaluated after truncation at either £ rcc = 10 or 
30. In both cases, following lEfstathioi] 120031) the input map was 
smoothed by a Gaussian beam with FWHM=7°. Smoothing im- 
proves the accuracy of the direct inversion method reconstructions, 
because it damps higher order modes. The input and output (recon- 
structed) coefficients a (c for the multipoles I = 2, . . . , 5 are com- 
pared in Fig. J3J and J4}. Only those modes with m = £ are shown 



because, as will be shown in Section llOl they (and the correspond- 
ing modes with m = -I) have the largest reconstruction uncertain- 
ties as a result of the removal of a large fraction of their power 
by the cut. Thus, they are the most sensitive to the method used 
in their estimation. Nevertheless, comparison of the other modes 
gives qualitatively similar results. 

Fig. 01 shows that the PE filter recovers the quadrupole and 
octopole amplitudes with a similar accuracy to the direct inversion 
method. However, for higher order modes, the PE filter performs 
better. This does not appear to be connected to the use of £ mc = 10 
for the direct inversion method; increasing this value to f rec = 30 
yields even worse reconstructions (Fig. J3)). The higher accuracy 
of the PE filter reconstructions is a consequence of explicitly ac- 
counting for the aliasing of higher order multipoles into the lower 
multipole range in the filter construction. 



7 APPLICATION TO THE WMAP DATA 

The PE filter was applied to WMAP maps with £ max = 30 and 
Aec = 10. We are mostly interested in the analysis of the lowest 
multipoles I = 2, . . . , 5 because these are the most sensitive to the 
cut and foreground emission. The monopole and dipole were not 
restored because they do not provide any relevant cosmological in- 
formation. In order to construct the filter matrix, the WMAP best- 
fitting ACDM model (fl A = 0.73, Q,„ = 0.27) with a running pri- 
mordial spectral was used. This choice is based on the fact that the 
spectrum provides a better fit to the WMAP data on large angular 
scales compared to a simple power-law spectrum. However, tests 
with a power-law spectrum yields similar results. The noise covari- 
ance matrices were constructed in the usual fashion, with pixel vari- 
ances determined by the number of observations N h s (p) in pixel p 
as provided by the WMAP team. 

The power spectrum will provide a useful complement to the 
studies of the multipole maps. We will use an estimator of the 
power spectrum given by 



c t = 



(16) 



2i+l 

where a (m denotes the restored harmonic coefficients. The defini- 
tion of the PE filter ensures that the estimator is unbiased. Given 
the error covariance matrix M, one can estimate the uncertainty in 
the quantity C f by 



{2i + iy 



^ M u fai&j - My) , 



(17) 



where the indices i,j in ^-ordering are expressed by 
i = £ 2 + £ + m + 1 . This estimation of the uncertainty takes 
into account only errors induced by the sky cut and noise; it does 
not take into account cosmic variance related errors which enter 
through the assumption of a cosmological power spectrum to build 
the PE filter. Table Q indicates the reconstruction accuracy of the 
method. 



7.1 The multipoles maps 

Maps of the recovered multipoles computed from the WMAP cor- 
rected co-added map are compared to the multipoles of the cleaned 
TOH and ILC maps in Fig. Because the details of the noise 
model are essentially negligible for the estimation of low-order 
multipoles, we have also used the PE filter as computed for the 
corrected co-added map to analyse the TOH and ILC maps. The 
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Table 1. Rms difference between the true and estimated power spectra 
^error ^ m ^^-2^ an( j me corre i a ti on coefficients ree of the f-order, m = ( 
modes estimated from a set of 10 000 simulations of the WMAP data anal- 
ysed after application of the Kp2 mask. In the third column, the signal to 
eiTor ratio is shown where Cc is the estimated power spectrum of the WMAP 
multipoles maps given by PE filtering. 



t 




C f /C«» 


r U 


2 


69 


2.8 


0.99 


3 


41 


12 


0.98 


4 


26 


9.6 


0.97 


5 


19 


17 


0.96 


6 


14 


5.2 


0.94 


7 


12 


11 


0.92 


8 


9.0 


6.3 


0.90 


9 


7.6 


7.6 


0.87 


to 


6.2 


7.3 


0.84 



comparison confirms the previous results of estimations of the mul- 
tipole amplitudes. Small discrepancies are seen only at I = 2, where 
the hot spots are closer to the Galactic equator and the amplitude is 
smaller for our results. The power spectrum of the maps is shown 
in Fig. J6j and summ arised in Table The cut sky analysis of 
Hinsh aw et alj 120031 shows subst antially less powe r at £ = 3,5 
and 7 than for the other estimates. Efstathiou (2003) has claimed 
that this is a consequence of the quadratic nature of this estima- 
tor. Our linear estimator results, and in particular the ILC results in 
Table J2j, are in good agreement with those of Efstathiou 1 2003). 

7.2 Dependence on the assumed angular power spectrum 

Because construction of the filters requires some assumption about 
the CMB anisotropy power spectrum CV, one should investigate 
how sensitive the amplitudes of the restored multipoles are to this 
choice and thus quantify the uncertainties in the amplitudes related 
to the assumption. To answer these questions we have applied fil- 
ters constructed using two other power spectra: the canonical one 
described above but with the square root of the cosmic variance 
either added or subtracted. In Fig. {7), we see that the power spec- 
trum of the resultant maps is changed only weakly by these pertur- 
bations, up to few percent. For the WMAP analysis presented here, 
the PE filtering method gives robust estimates of the amplitudes of 
the low-order multipoles, essentially independent of the assumed 
input power spectrum. 

7.3 Dependence on frequency 

Comparison of the multipoles maps determined from the fore- 
ground corrected WMAP data frequency-by-frequency allows us to 
determine the extent to which the derived structures are cosmologi- 
cal and the degree of residual foreground contamination. In Fig. {S}, 
it is seen that the foreground cleaning of the A"-band map gives the 
worst results. The multipole maps differ substantially from the rest 
of the maps showing excess of the power near the Galactic plane. 
Nevertheless, apart from the quadrupole, there is a qualitative sim- 
ilarity in the distribution of hot and cold spots relative to the other 
frequency channels. In particular, the octopole is planar with ap- 
proximately the same preferred axis, along which the power is sup- 
pressed, for all frequencies. The quadrupole is also planar, although 
changes of the plane are more pronounced, especially for the K and 
Ka bands. Small deviations between the Q, V and W bands indi- 
cates that, after foreground cleaning, these maps are dominated by 



the CMB anisotropy component and can be treated as a reliable 
picture of the CMB multipoles. 

The corresponding power spectra shown in Fig. and sum- 
marised in Table J3J clearly show that the applied foreground cor- 
rection does not remove all the contamination for the K and (par- 
tially) Ka bands. This suggests that the Kp2 cut is not sufficient to 
remove the strong signals near to the Galactic plane at those fre- 
quencies (despite the fact that the Kp2 mask is derived from con- 
siderations of the thresholded temperature distribution of the K- 
band map). The strong synchrotron and free-free emission extends 
to higher Galactic latitudes at these lower frequencies. 

7.4 Dependence on the sky cut 

More substantial changes of the multipoles are seen with then ap- 
plication of more aggressive sky cuts. In Fig. <10> . we show the 
multipoles determined from the WMAP corrected co-added data af- 
ter masking by the Kp2, KpO, 20+ and 30 masks. All multipoles 
lose their characteristic features for mask 30. Nevertheless, up to 
sky cut 20+ the multipoles £ = 3 and t = 5 do not change signif- 
icantly. The multipoles £ = 2 and I = 4 are more sensitive to the 
mask and change substantially even for smaller cuts. 

The corresponding power spectra are presented in Fig. jl H 
and Table It should be noticed that generally the uncertainties 
increase rapidly with the width of the sky cuts. Because our esti- 
mation of the errors J 171 depends on the multipole power, the few 
exceptions from this rule, £ = 6 for the 20+ mask, £ = 7, 8, 9 for the 
30 mask, are a consequence of the substantially lower amplitudes 
of the multipoles estimated in those cases. Notice that for very low 
power, the estimation of the square of the error can be negative. 
Because we have such a situation for the quadrupole estimated for 
the cut 30, we do not quote the corresponding uncertainties in the 
table. 

The general trend for increasing uncertainties with increasing 
sky cut is easily seen in Fig. <12> where linear regressions, anal- 
ogous to that in Section [5] is shown for various masks. It should 
be noticed that uncertainties in the estimation of the multipoles 
i = 4, 5 is especially large for the cut 30. The true and reconstructed 
coefficients are nearly uncorrelated. Slightly better correlations are 
seen for modes with m 4- —£, £, though the trend is still maintained. 
The poor correlation between input and recovered coefficients im- 
plies that for I > 5 their estimation is no longer credible. Indeed, 
the large uncertainties associated with the widest cut 30 helps to ex- 
plain why the appearance of the estimated WMAP multipoles is dif- 
ferent from that estimated with narrower masks. Nevertheless, the 
multipole amplitudes remain statistically consistent over the range 
of masks at approximately the 2<x level. 



8 UNCERTAINTIES OF THE MULTIPOLES MAPS 

One of the merits of the PE filtering method is the possibility to 
directly estimate the uncertainties in the derived results, and the 
preservation of the well-defined noise properties of the original 
maps. This is a clear advantage of the method over the MEM ap- 
plied by the WMAP team I Ben nett et alj|2003al which results in 
anisotropy maps with computationally complex noise properties 
which are largely inappropriate for most CMB analyses. Indeed, 
it is generally true that the noise evaluation for linear methods are 
computationally more feasible than for non-linear methods such as 
MEM. Here we derive errors both in pixel and spherical harmonic 
space. 
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Table 2. Power spectrum amplitudes AT 2 = (if + l)Cc/2n (in fiK 2 ) for the foreground corrected WMAP frequency maps and related combinations estimated 
by application of the PE filter for the Kp2 sky coverage. Q+V+W denotes the corrected co-added map and * denotes estimation on the full sky. The direct 
inversion method was applied to the corrected co-added map for the Kp2 mask. Uncertainties do not take into account cosmic variance. 



AT; AT; AT; AT; AT; AT; AT; AT; AT;. 

2 3 4 5 6 7 5 9 l( 



K 963 + 79 1710+ 160 3290±220 2070+ 130 710+ 110 1060 + 120 1320+ 120 1330+ 150 1680+ 180 

Ka 176 + 26 1060+ 120 1190+ 110 1550+ 110 312 + 55 960+ 100 786 + 84 920+ 130 1100+ 140 

Q 158 + 32 990+ 110 875 + 80 1610+ 120 459 + 77 1100 + 110 683 + 75 880+ 130 840+ 120 

V 206 + 40 930+ 100 755 + 63 1600+ 120 523 + 79 1070 + 110 644 + 69 800+ 130 750+ 110 

W 219 + 42 910+ 100 804 + 67 1580+ 130 553 + 79 1110 + 120 645 + 69 780+ 130 770+ 110 

Q+V+W 182 + 36 950+ 110 804 + 70 1600+ 120 491 + 78 1090 + 110 656 + 71 830+ 130 790+ 110 

Direct inversion 322 + 36 873 + 96 1530+ 120 1370+ 110 569 + 78 870+ 100 1037 + 99 950+ 140 1090+ 130 

1LC 248 + 45 950+ 110 750 + 61 1620+ 130 615 + 85 1120 + 120 619 + 66 790+ 130 639 + 98 

TOH 279 + 49 870 + 98 613 + 49 1420+ 110 586+ 89 1040 + 110 631 +63 860+ 140 760+ 110 

ILC* 195 +40 1050+ 110 834 + 66 1670+ 130 606+ 81 1290 + 140 660 + 68 700+ 130 780+ 110 

TOH* 202 + 41 866 + 99 651 + 50 1350+ 110 608 + 90 1010 + 110 681 +72 710+ 120 890+ 120 



Table 3. Power spectra of the multipoles AT 2 = ((£. + \)Ccl2n (in fiK 2 ) 
estimated by the PE filtering after applying various masks. Uncertainties do 
not take into account cosmic variance. 







Mask 






Kp2 


KpO 


20+ 


30 


1 = 2 


182 ± 36 


294 ± 88 


200 ± 130 


77 


f = 3 


950 ± 110 


870 ± 170 


890 ±310 


830 ± 380 


e = 4 


804 ± 70 


723 ± 94 


970 ± 280 


790 ±310 


£ = 5 


1600 ± 120 


1470 ± 170 


820 ± 280 


980 ± 450 


AT 2 € = 6 


491 ± 78 


480 ± 110 


200 ± 54 


690 ± 260 


1 = 1 


1090 ± 110 


1110 ± 190 


1460 ± 380 


700 ± 240 


f =8 


656 ±71 


750 ± 130 


920 ± 250 


560 ± 130 


e = 9 


830 ± 130 


780 ± 180 


860 ± 250 


750 ± 220 


c=w 


790 ± 110 


610 ± 140 


590 ± 140 


730 ± 170 



8.1 The maps of the standard deviation 

To study the spatial distribution of the uncertainties of the restored 
multipole maps, we use the simulations of the WMAP corrected 
co-added map as described in Section|4] The map 

A 2 (p)e rrol - = 2 Z (CC) (P)^»M - (18) 

I i.ni! f 2 .'»2 

shows the spatial distribution of the squared standard deviation of 
the difference between the restored and true multipole maps. fl" r ™ 
are the corresponding decomposition coefficients of the difference 
maps and terms (flf"™^™^} are averaged over 10 000 simulations. 
Using the identity 



Jf 1 ,m 1 (p)^ 2> m 2 (P) 



I 

(.III 



(2d + l)(2£ 2 + l)(2€ + 1) 



An 



xi ,„. ,„ )( ,» i>-""- 

where is the Wigner 3i symbol, the expression 

1 m\ m 2 m ) a j j t- 



i 1 Si can be written as 



A 2 ( P )e, 



^ b C.mY(,m(P) . 



where b tm are the appropriate coefficients. The A 2 rrOT map is 
composed of modes < £ < 2t\, because £ obeys the triangular 
condition {l = t\ +t 2 Ji +4-1,..., \( x - £ 2 \) and t x = t 2 . 



The A 2 . ror maps computed for the Kp2 mask are shown in Fig. 
J 1 3i in logarithmic scale. As expected, as a result of the lack of 
the data in the cut part of the sky the reconstructed multipoles have 
the largest uncertainties in this region, especially close to the Orion 
and Ophiuchus exclusion zones. The multipoles C = 2 and 3 have 
the smallest uncertainties inside the cut. This is a consequence of 
the fact that the multipole maps characterize the CMB anisotropy 
on angular scales ~ If this is bigger than the angular width, 
A9, of the mask then the excised region of the multipole map can 
be easily extrapolated from the data outside of the cut. Otherwise, 
there will be fluctuations with typical angular size smaller than A9, 
which will be completely contained within the cut and their recon- 
struction from the data will be much less accurate. The statistical 
isotropy assumption used in our approach is then very helpful to 
reduce uncertainties in these cases. For the Kp2 mask, the widest 
part, cutout toward the Galactic centre, has width Ad ss 40° . 

The A 2 rror map corresponding to the sum of the multipoles 
i = 2, ... ,5 shows much bigger uncertainties than for the corre- 
sponding individual multipoles. This is a result of the additional 
contributions from the correlations between errors in the multi- 
poles. These correlations are easily seen in the covariance matrix 
M. 



8.2 The covariance matrix 

The map of the standard deviation does not show all of the infor- 
mation about the errors. The covariance matrix M shows that the 
uncertainties are correlated in the spherical harmonic space. 

As mentioned in Section|2| the covariance matrix, M, consists 
of three components: the noise component and, induced purely by 
the cut, contributions from correlations between the restored mul- 
tipoles and correlations with higher order multipoles. In Fig. J I -41 
we show that part of the covariance matrix corresponding to the 
multipoles £ = 2, . . . , 5 for the PE filtering of the WMAP corrected 
co-added map, together with the corresponding sub-matrices for 
the separate components. The largest contribution arises from cor- 
relations with higher order modes. The correlations between the 
low-order reconstructed multipoles are weaker although still non- 
negligible. The smallest contribution comes from the noise com- 
ponent, which is plotted multiplied by a factor of 50 to render it 
visible. The approximately axisymmetric cut causes stronger cor- 
relations between modes with £ value different by 2. Notice that the 
modes with the biggest reconstruction uncertainties are for m = —t 
and m = £, which is a consequence of removing a large fraction 
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of the power of these modes, mainly concentrated in the Galactic 
plane, by the cut. 



9 FOREGROUND SEPARATION 

In this section, we will reconsider the sensitivity of the estimated 
CMB multipoles to the Galactic foregrounds. To achieve this, we 
derive a CMB anisotropy map by the spectral decomposition of 
multi-frequency WMAP data by making assumptions about the fre- 
quency dependence of the foreground components. The advantage 
of this m ethod is that, unlike with the template fits o flBennett et alj 
(2003b), we assume nothing about the morphological distribution 
of the foreground emission. Thus foregrounds which are not well 
traced by the existing temp lates, but with sim ilar spectra, e.g. the 
free-free haze proposed by iFinkbeiner (2003), are taken into ac- 
count. The disadvantage is that, to separate the components prop- 
erly, we need to know their spectral dependence. Given that our 
intention is not to provide a definitive foreground separation anal- 
ysis, we have considered a broad range of spectral index values 
to aggressively test the dependence of the reconstructed low-order 
multipoles on any foreground assumptions. Nevertheless, it should 
be recognised that our assumptions may lead to worse foreground 
subtraction than in the WMAP analysis, although this is unlikely 
given that the range of spectral properties adopted here are well 
motivated by previous results in the literatur e. 

We consid er a method followin g jBennet t et alj ll992h and the 
formalism from Dodelson & Stebbins 1 1994). We suppose now that 
we observe the sky in m different frequency channels, so that at 
each point r on the sky we measure m different temperatures d,(r), 
i = l,...,m. Let us assume also that we know the temperature 

spectra fj(v) of all n components Xj(r), j = 1 n. Then we can 

write 

d(r) = F • x(r) + n(r) , (19) 

where the vector n(r) corresponds to the instrumental noise in the 
various channels and F (not to be confused with the filter matrix 
used in Section[2) is the mxn frequency response matrix given by 

Fij= f w,(v)fi(v)dv , (20) 
Jo 

w,-(v) being the frequency response (bandwidth) of the ith channel. 
Here, we will approximate w,(v) by Dirac's delta w,(v) = S(v - v,) 
where v, corresponds to the mean frequency in the ith channel, so 
that Fjj = fj(vj). If the number of frequencies, m, equals the num- 
ber of components n and noise is negligible, the separation of the 
components can be obtained by inverting the frequency response 
matrix, 

x(r) = F~' • d(r). (21) 

For the CMB component, if temperature is expressed in thermody- 
namic units, F;i = 1 over all frequencies as a result of the black- 
body nature of the spectrum. For the range of frequencies of interest 
here, the foreground spectra may be reasonably well modelled by 
a power law in antenna temperature, fj(v) oc iA , where the spec- 
tral index y8 ; depends on the component. To express the spectral 
dependence in thermodynamic units, the function fj(v) must in- 
clude the frequency dependent antenna to thermodynamic temper- 
ature conversion factor, p(v). Thus, finally the functions take the 

form fj(y) = p(v)(v/v jj ' , where v 0j is the reference frequency 
for component j at which the foreground amplitude is normalised. 



For the WMAP data, sky maps are available at 5 frequen- 
cies: K, Ka, Q, V, and W bands 7 (i = 1,...,5, respectively) 
and these may be interpreted as consisting of 5 components: the 
CMB anisotropy, Galactic thermal dust, synchrotron, free-free and 
anomalous dust emissions (j = 1, . . . , 5, respectively). 

The synchrotron, free-free and anomalous dust components 
will be searched at A"-band frequency v 0s = Vgjf = v 0cl = 
23 GHz and the thermal dust component at W-band frequency 
vo.d = 94 GHz. It should be recognized that the Galactic emis- 
sion does in general have spectral properties that depend on po- 
sition on the sky (especially in the the Galactic plane), however 
we will neglect these variations and use a fixed index over the 
high-latitude parts of the sky we are interested in. This is no dif- 
ferent to the assumption that the template fits can be described by 
single scale factors over the same region. Even in the case of the 
WMAP ILC method, the high-latitude emission (as defined by the 
Kp2 mask) was treated as a single contiguous region over which the 
foregrounds were assumed to exhibit unvarying spectral behaviour. 
Only within the Galactic plane as defined by the Kp2 mask were 
additional regions defined to allow for spectral variations. 

For the thermal dust compon ent the spectral ind ex generally 
lies in the range 1.6 < fid < 2.5 (Dupac et al. 2002). Th r ee val - 
ues were selected: /3 rf = 1.7 suggested by Finkbeiner et al. 1 1999), 
Pj = 2.2 which is the mean value of the index determined from 
the MEM-fitting dust component derived from the WMAP data by 
iBennett et alj l2003bl) . and the upper limit of the range fi d = 2.5. 

The synchrotron spectral index is negative and again three val- 
ues were selected: f} s = -3.3 which approximately corresponds to 
the mean spectral index between the Ka and Q bands of the WMAP 
MEM-fitting synchrotron component, p s = -3.1 which is a better 
fit to the Galactic halo emission, and f} s = -2.9 more appropriate 
for the Galactic plane. 

The spectral depen dence of the so -called anomalous dust com- 
ponent was studied by iBandav et alj 120031) . Over the frequency 
range 19-90 GHz, the anomalous dust has a thermal dust-like mor- 
phology but a synchrotron-like spectrum with an index /3„ » -2.5. 
Unfortunately, the similarity of the latter two spectra results in the 
the frequency response matrix becoming poorly conditioned and 
the clean separation of the components is not possible. To sidestep 
thi s issue, we ado pted the spinning dust model spectrum proposed 
by Draine & Lazarian 1 1998) (hereafter DL98) 8 . In particular, we 
consider a three-component model of the spinning dust containing 
contributions from dust found in the cold neutral medium (CNM), 
warm neutral medium (WNM) or warm ionised medium (WIM). 
The preferred DL98 model has these components in the ratio 0.43, 
0.43 and 0.14, respectively. It should be noted that the model pre- 
dictions are fundamentally uncertain because changes in the grain 
composition can affect the turn-over frequency, the spectral nor- 
malization, and so forth. In addition, variations in the mixture of 
contributions from the three phases of the interstellar medium can 
affect the integrated spectrum, making our fit a very poor approxi- 
mation. We have tested the sensitivity of the final CMB multipole 
reconstruction to the assumptions by replacing the preferred model 
mix by either the CNM, WNM or WIM spectra, respectively. The 
results are remarkably stable to these variations, the angle between 
the quadrupole and octopole axes changing by less than 3°. Fur- 



7 Note that, before proceeding with the separation, these were smoothed to 
a common angular resolution of 1 ° . 

8 Detailed spectra are available at 

http://www. astro. princeton. eduj draine/dust/dust. mwave. spin.html. 
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thermore. Ide Oliveira-Costa et alJ J2004b) find spectral behaviour 
in data from 10-90 GHz in reasonable agreement with the predic- 
tions of DL98, and we persist with these assumptions. 

Nevertheless, further degeneracy in the components caused by 
the similarity of the free-free spectral p/f as -2.15 and the syn- 
chrotron indices, renders it impossible to perform a robust sepa- 
ration of all 5 components. We have to eliminate the free-free or 
synchrotron component from the linear system. Because the for- 
mer is better known (in the sense that there are less variations in the 
spectral index over the sky) and of lesser contribution to the fore- 
ground signal, we adopt a free-free template to describe its spatial 
morphology and subtract it from the 5 temperature maps, d(r), with 
the correlation coefficients derived by Ben nett et al .(2003b) (as de- 
scribed in Section[3J. However, a major uncertainty when using the 
Ha template is the absorption of the Ha by foreground dust. To 
study the sensitivity of the lowest multipoles to the dust absorption 
correction we wil l use two template s of the free-free emission: Ho- 
template given bv lFinkbeine3 f2003) (it was used by iBennett et alJ 
2003b) uncorrected for dust absorption and Ha template corrected 
for dust absorption given by iDickinson. Davies & Davisl 120031) . 
Unfortunately, dust may not be th e only problem in removing free- 
free emission. Finkbeiner (2004) has suggested that the Ha map 
does not trace all of the free-free emission, and proposed the ex- 
istence of a so-called free-free haze which has the same spectral 
dependence as free-free component. The influence of the free-free 
haze on the estimation of the lowest multipoles will be studied us- 
ing maps corrected and uncorrected for this component. 

After elimination of the free-free component, an overdeter- 
mined linear system of 5 equations and 4 unknowns remains to be 
solved. This is performed using the singular value decomposition 
(SVD) method, because it gives the minimum variance solution. 
Then, the PE filter was applied to the derived CMB component 
maps. Residuals associated with poorly subtracted Galactic plane 
emission were circumvented by application of the Kp2 mask. 

The influence of the dust absorption and free-free haze cor- 
rections on the lowest WMAP multipoles is shown in Fig. d 1 5i . The 
dust and synchrotron spectral indices were fixed here as f3j = 2.2 
and f} s = -3.1, respectively. One can see that quadrupole is espe- 
cially sensitive to these corrections. The dust absorption correction 
shifts the hotter spots of the quadrupole toward the Galactic poles 
and the free-free haze correction shifts the spots toward the Galactic 
plane. Curiously, when both corrections are applied, the quadrupole 
remains largely unchanged. Notable changes are seen also for the 
£ = 4 multipole, but the results for £ = 3 and 5 are very stable with 
respect to the corrections. 

The dependence of the 1 = 2 multipole on the thermal dust and 
synchrotron spectral indices is shown in Fig. < Kit . We do not show 
analogous figures for I = 3,4,5 because they remain unperturbed 
it = 3, 5) or show minimal changes (I = 4) with respect to the spec- 
tral indices. It should be noted that an increase in the synchrotron 
spectral index tends to shift the hot spots of the quadrupole toward 
the Galactic poles and enhance its amplitude. Conversely, increas- 
ing the dust spectral index shifts the spots toward the Galactic plane 
and decreases their amplitude. Thus, for the flattest dust fid = 1.7 
and synchrotron fi s = -2.9 indices, the quadrupole hot spots are at 
the Galactic poles. Similar changes though much smaller are also 
seen for multipole t = 4. 

The stability of the multipoles £ = 3 and 1 = 5 with respect 
to the spectral indices, combined with the earlier results related to 
the removal of the free-free template, is compelling evidence in 
favour of the claim that the planarity of the octopole is not caused 
by the Galactic foregrounds. Conversely, the strong dependence of 



Table 4. Galactic coordinates of the quadrupole and octopole axes for dif- 
ferent masks. 









Kp2 


KpO 


20+ 


e = 


2 


Longitude 


273° ± 8° 


298° ± 5° 


276° ± 53° 






Latitude 


68° ± 3° 


40° ± 16° 


34° ± 28° 


e = 


3 


Longitude 


307° ± 4° 


301° ± 6° 


267° ± 22° 






Latitude 


63° ± 2° 


62° ± 3° 


66° ± 9° 


Am 


>le between axes 


15° ± 3° 


22° ± 14° 


33° ± 23° 



Table 5. Galactic coordinates of the quadrupole and octopole axes for dif- 
ferent frequencies. In all cases the multipoles were derived from input maps 
where the Kp2 mask was applied. The coordinates have uncertainties as in 
Table |4j. 





K 


Ka 


Q 


V 


W 


1 = 2 Longitude 


107° 


228° 


262° 


282° 


280° 


Latitude 


1° 


29° 


66° 


63° 


71° 


t = 3 Longitude 


299° 


312° 


310° 


307° 


306° 


Latitude 


74° 


69° 


65° 


63° 


62° 


Angle between axes 


105° 


61° 


19° 


11° 


13° 



the quadrupole amplitude and morphology on the spectral prop- 
erties of the Galactic foreground emission implies that the latter 
must be known in some detail in order to recover the quadrupole 
anisotropy with some degree of precision. Indeed, this will also 
have important consequences for the statistical significance of its 
alignment with octopole. This issue is studied in next section. 



10 THE QUADRUPOLE AND OCTOPOLE PREFERRED 
AXES 

As initially noticed by Ide Oliveira-Costa et alJ l2004al) . and seen 
here in Fig. the quadrupole and octopole are rather planar, with 
most of their hot and cold spots placed on a single plane in the sky. 
de Oliveira-Costa et al. (2004a) use a wave function formalism to 
quantify a preferred axis for the multipoles and study the signifi- 
cance of their alignment. However, their analysis was based on the 
full sky TOH map which may contain residual contamination in the 
Galactic plane region. We repeat the analysis here on the multipoles 
derived from the cut sky analysis (thus avoiding this problem), and 
quantify the additional statistical uncertainties introduced by the 
mask. 

Given the covariance matrix M, we can estimate uncertain- 
ties of the preferred axes coordinates on the sky. For each of the 
masks: Kp2, KpO and 20+ the covariance matrix M was used to per- 
form 100 constrained realizations of the WMAP quadru pole and oc- 
topole. Then, using the same technique as ide Oliveira-Costa et alJ 
(2004a) the preferred axis for each of the multipoles was found. 
The uncertainties in the axis coordinates are summarised in Table 
J4j. Notice that the multipoles axes coordinates substantially vary 
for different masks. The most changeable is the quadrupole, partic- 
ularly in latitude, whereas the octopole axis coordinates are much 
more robust. 

Table J5} indicates the frequency dependence of the preferred 
axis for the quadrupole and octopole. One should notice that the 
best foreground corrected octopole maps (those of the Q, V and W 
bands) are stable and prefer approximately the same axis: (I, b) ~ 
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Table 6. Galactic coordinates of the quadrupole and octopole axes for four 
variants of the free-free emission correction of the WMAP maps (see text). 
Columns (1), (2), (3) and (4) correspond to, respectively: no free-free haze 
correction and no dust absorption correction, free-free haze correction and 
no dust absorption correction, no free-free haze correction and dust absorp- 
tion correction, free-free haze correction and dust absorption correction. 
The dust and synchrotron spectral indices are fi c [ = 2.2 and p s = —3.1, 
respectively. 





No dust 


corr. 


Dust corr. 




(1) 


(2) 


(3) (4) 


t — 1 Longitude 


297° 


292° 


300° 296° 


Latitude 


53° 


69° 


33° 62° 


C = 3 Longitude 


306° 


301° 


301° 299° 


Latitude 


62° 


62° 


62° 60° 


Angle between axes 


10° 


8° 


29° 2° 



Table 7. Galactic coordinates of the quadrupole axes for various spectral 
indices of the dust and synchrotron emission. The first number in brackets 
is longitude and the second latitude. In all cases the maps were masked by 
Kp2 mask, thus the coordinates have uncertainties the same as in Table J4} 
for Kp2 mask. 



/^synch 




-3.3 


-3.1 


-2.9 


1.7 

Ai ust 2.2 
2.5 


(298°, 59°) 
(291°, 68°) 
(287°, 71°) 


(298°, 45°) 
(296°, 62°) 
(294°, 65°) 


(297°, 25°) 
(298°, 43°) 
(297°, 53°) 



11 APPLICATION TO THE COBE-DMR DATA 

As a useful cross-check, in particular against systematic effects, we 
investigate whether the lowest order multipole maps derived from 
COBE-DMR are consistent with the WMAP results. We treated 
the COBE-DMR data in similar way to the WMAP data. Firstly, 
we corrected the maps for the Galactic foreground using the 20+ 
mask and the Gala ctic emission templates as was described in 
iGorski et alj ll99cj) and then co-added the 53 and 90 GHz maps 
using inverse-noise- variance weights. The 31.5 GHz data were not 
used because under such a weighting scheme the contribution of 
these least sensitive maps is minimal. The PE filter was constructed 
using the best fit to the data power spectrum p arameterized by 
bond & EfstathioJl987llFabbri. Lucchin & Matarresel 19871) 

r _/ 2 \_ n 2 4n T(e + ^)r( 9 -f) 

U < - \ a C,ml - Unm-PS g ^ + 5^1-^3+^ ' ( - ZZJ 

where T denotes the gamma function. The rms quadrupole normal- 
isation Q rms -ps and spectral index n were estimated in Gors ki et alJ 
Jl99fj) . Qrms-ps = 15.3+j-g fiK, n= 1.2 ± 0.3, and we used the val- 
ues of the parameters which maximize the likelihood function The 
noise covariance matrix was estimated analogously to the WMAP 
matrix. i rec and C mlix were the same as for the WMAP PE filter. 

The COBE-DMR and WMAP multipoles maps for the same 
20+ mask are compared in Fig. 1171 . One sees that the maps are 
similar. The COBE-DMR octopole appears to be planar like the 
WMAP octopole and prefers roughly the same axis, in the direc- 
tion of (l,b) ~ (277°, 63°). The hot spots of the COBE-DMR 
quadrupole are shifted more toward the Galactic poles than the 
spots of the WMAP quadrupole, which may reflect differences in 
the foreground removal. 



(308°, 63°) with a deviation of order 2°. The quadrupole axis is 
more poorly determined. Significant differences between the Q, V, 
W and K, Ka frequency bands suggest the presence of residual fore- 
grounds at the lowest frequencies. 

As was shown in Section [5] the quadrupole of the CMB 
map derived by the linear projection depends significantly on tem- 
plates used to eliminate the free-free emission and assumptions 
concerning the spectral dependence of the dust and synchrotron 
emission. The coordinates of the quadrupole and octopole pre- 
ferred axes for variants of the free-free model emission are shown 
in Table ^6). Again, the octopole is much more stable than the 
quadrupole. Correction for the putative free-free haze emission 
shifts the quadrupole axis toward the Galactic poles (larger lati- 
tude), and the dust absorption correction toward the Galactic plane 
(smaller latitude). 

The octopole axis is also stable and does not change for dif- 
ferent spectral indices of the dust and synchrotron emission. For all 
studied indices it is (/, b) = (299°, 60°). The quadrupole axis, as in 
previous cases, is much less stable (see Table 0). 

In summary of this section it is worth emphasizing that the 
above results show a much better determination of the octopole and 
its preferred axis than for the quadrupole. Significant uncertainties 
of estimation for the quadrupole preferred axis and its strong de- 
pendence on the details of templates used to eliminate free-free 
emission and foreground spectral indices must increase the low 
probability of the WMAP quadrupole-octopole alignment estimated 
bv lde Oliveira-Costa et alJ l2004al) . The importance of quantifying 
errors as a result of the foreground subtraction should be clear. 



12 SUMMARY 

In this paper, we have re-examined the nature of the low-order 
multipoles derived from the WMAP data. Previous work has sug- 
gested that these modes demonstrate unusual properties, but such 
results have been derived from full sky analyses of maps of the 
CMB anisotropy derived by methods which are susceptible to resid- 
ual Galactic foregrounds (Eriksen et al. 2004a), particularly in the 
Galactic plane. By application of a PE filter on the cut sky we cir- 
cumvent these difficulties, and have estimated the amplitudes of 
the £ = 2 — 5 multipoles and their dependence on frequency, sky 
cut and foreground correction. Because the construction of the fil- 
ter takes into account correlations of the lower order modes with 
higher order modes, it significantly improves the accuracy of mul- 
tipole e stimation in com parison with the direct inversion method 
used bv lEfstathioul <2003l) . 

Our studies of the dependence on frequency of the multipole 
amplitudes showed that: 

(i) the Kp2 mask is insufficient to remove strong foreground con- 
tamination near the Galactic plane for the K and Ka frequency 
bands even for the foreground corrected WMAP data; 

(ii) small deviations between the I = 3, 4, 5 multipole deter- 
mined from the Q, V and W frequency bands indicates that, after 
foreground cleaning, these multipoles are dominated by the CMB 
anisotropy component outside the Kp2 mask; 

(iii) the octopole preferred axis, approximately (/, b) ~ (308°, 63°), 
shows deviations of order two degrees for the Q, V and W fre- 
quency bands; 
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(iv) noticeable changes of the quadrupole suggest a non-negligible 
residual foreground emission. 

That the quadrupole (and to a lesser extent I = 4) may be 
compromised by such residuals is supported by the changes ob- 
served after the application of increasing sky cuts. It also appears 
that the estimation of the quadrupole is sensitive to the dust ab- 
sorption correction of the free-free template and free-free haze cor- 
rection. This may provide yet further independent confirmation of 
significant contamination of the quadrupole by residual Galactic 
foregrounds. 

Further tests of the sensitivity of the multipoles to the spec- 
tral dependence of the foreground emission have utilised a linear 
projection technique for foreground removal. Studies of the depen- 
dence of the multipole amplitudes and morphologies on the dust 
and synchrotron spectral indices showed that: 

(i) multipoles C = 3,4,5 are independent of the indices and do 
not vary with respect to them; 

(ii) increase of the synchrotron spectral index shifts hotter 
spots toward the Galactic poles and enhances amplitude of the 
quadrupole; 

(iii) increase of the dust spectral index shifts hotter spots toward 
the Galactic plane and decreases amplitude of the quadrupole. 

In conclusion, it remains premature to make definitive state- 
ments in relation to the quadrupole. Large uncertainties of the 
quadrupole preferred axis should substantially increase the very 
low probabilities of the WMAP quadrupole-octopole alignment es- 
timated by de Oliveira-Costa et al. (2004a) and other quadrupole 
correlations st udied bv ISchwarz et alJ 120041) . To similar conclu- 
sion came also lSlosar & Seliald 12004) who showed that much of 
the alignment between quadrupole and octopole is the result of 
foreground contamination. Nevertheless, estimates of £ = 3,4,5 
are quite robust, and consistent with previous values. The observed 
anisotropy structures are also present for in COBE-DMR data, and 
are therefore unlikely to be associated with instrumental systematic 
artifacts. 



13 ACKNOWLEDGMENTS 

We a cknowledge the use of CMBFAST fceliak & Zaldarriagj 
1996). So me of the results in th is paper have been derived using the 
HEALPix I Gors ki et alJll998f) software and analysis package. We 
acknowledge the use of the Legacy Archive for Microwave Back- 
ground Data Analysis (LAMBDA). Support for LAMBDA is pro- 
vided by the NASA Office of Space Science. PB acknowledges the 
financial support provided through the European Community's Hu- 
man Potential Programme under contract HPRN-CT-2000-00124, 
(CMBNET), and Max-Planck-Institut fuer Astrophysik in Garch- 
ing where parts of this work were carried out. 



REFERENCES 

Banday A.J., Dickinson C, Davies R.D., Davis R.J., 2003, MN- 
RAS, 345, 897 

Bennett C.L., Smoot G.F, Hinshaw G, et al. 1992, Apl, 396, L7 
Bennett C.L., Banday A.J., Gorski K.M., et al. 1996, ApJ, 464 , 
LI 

Bennett C.L., Halpern M., Hinshaw G, et al. 2003a, ApJS, 148, 1 
Bennett C.L., Hill R.S., Hinshaw G., et al. 2003b, ApIS, 148, 97 
Bond J.R., Efstathiou G. 1987, MNRAS, 226 , 655 



Bunn E.F, Fisher K.B., Hoffman Y., Lahav O, Silk I., Zaroubi S., 

1994, ApJ, 432, L75 
Copi C.J., Huterer D., Starkman G.D., 2004, Phys. Rev. D, 70, 

043515 

Cornish N.J., Spergel D.N., Starkman G.D., Komatsu E., 2004, 

Phys. Rev. Lett., 92, 201302 
de Oliveira-Costa A., Tegmark M., Zaldarriaga M., Hamilton A., 

2004a, Phys. Rev. D, 69, 063516 
de Oliveira-Costa A., Tegmark M., Davies R. D., Gutierrez CM., 

Lasenby A.N., Rebolo R., Watson R.A., 2004b, ApJ, 606, 89 
Dickinson C, Davies R.D., Davis R.J., 2003, MNRAS, 341,369 
Dodelson S„ Stebbins A., 1994, ApJ, 433, 440 
Draine B.T., Lazarian A., 1998, ApJ, 494, L19 
Dupac X., Giard M., Bernard J.-P, et al. 2002, A&A, 392, 691 
Efstathiou G., 2003, MNRAS, 346, L26 
Efstathiou G., 2004, MNRAS, 348, 885 

Eriksen H.K., Banday A.J., Gorski K.M., Lilje P.B., 2004a, ApJ, 
612, 633 

Eriksen H.K., Hansen F.K., Banday A.J., Gorski K.M., Lilje P.B., 

2004b, ApJ, 605, 14 
Fabbri R., Lucchin F, Matarrese S., 1987, ApJ, 315 , 1 
Finkbeiner D.P., 2003, ApJS, 146 , 407 
Finkbeiner D.P., 2004, ApJ, 614, 186 

Finkbeiner D.P., Davis M., Schlegel D.J., 1999, ApJ, 524 , 867 
Gould A., 1993, ApJ, 403 , L51 
Gorski K.M., 1994, ApJ, 430 , L85 

Gorski, K.M., 1997, in Bouchet, F.R., Gispert R., eds, Proceedings 
of the XVIth Moriond Astrophysics Meeting, Microwave Back- 
ground Anisotropies, Editions Frontieres, Gif-sur-Yvette, p. 77 
Gorski K.M., Banday A.J., Bennett C.L., Hinshaw G, Kogut A., 

Smoot G.F, Wright E.L., 1996, ApJ, 464 , Lll 
Gorski K.M., Hivon E., Wandelt B.D., 1998, in Banday A.J., Seth 
R., Da Costa L., eds, Proceedings of the MPA/ESO Cosmology 
Conference, Evolution of Large Scale Structure. PrintPartners 
Ipskamp, NL, p.37 
Haslam C.G.T., Klein U., Salter C.J., Stoffel H., Wilson W.E., 
Cleary M.N., Cooke D.J., Thomasson P., 1981, A&A, 100 , 209 
Hinshaw G, Spergel D.N., Verde L., et al. 2003, ApJS, 148, 135 
Kogut A., Banday A.J., Bennett C.L., Gorski K.M., Hinshaw G., 

Smoot G.F, Wright E.L., 1996, ApJ, 464 , L5 
Larson D.L., Wandelt B.D., 2004, ApJ, 613, L85 
Luminet J-P, Weeks J., Riazuelo A., Lehoucq R., Uzan J-Ph., 

2003, Nat, 425, 593 
Maartens R„ Ellis G.F.R., Stoeger W.R., 1996, A&A, 309 , L7 
Schwarz D.J., Starkman G.D., Huterer D., Copi C.J., 2004, 

preprint (astro-ph/0403353 I 
Seljak U., Zaldarriaga M., 1996, ApJ, 469, 437 
Slosar A., Seljak U., 2004, Phys. Rev. D, 70, 083002 
Smoot G.F, Bennett C.L., Kogut A., et al. 1992, ApJ, 396 , LI 
Spergel D. N., Verde, L., Peiris, H. V, et al. 2003, ApJS, 148, 175 
Stark Ph.B., 1993, ApJ, 408 , L73 

Tegmark M., de Oliveira-Costa A., Hamilton A.J., 2003, Phys. 

Rev. D, 68, 123523 
Vielva P., Martinez-Gonzalez E., Barreiro R.B., Sanz J.L., Cayon 

L., 2004, ApJ, 609, 22 
Wiener N., 1949, Extrapolation, interpolation and smoothing of 

stationary time series with engineering applications. The MIT 

Press, Cambridge, MA, USA 
Zaroubi S., Hoffman Y, Fisher K.B., Lahav O, 1995, ApJ, 449, 

446 



12 P. Bielewicz, KM. Gorski and A.J. Banday 



APPENDIX A: COVARIANCE MATRIX 

The posterior probability is a multivariate Gaussian of the form 

Pr(a £ |c) cc exp (-(a £ - a x ) r ■ M" 1 ■ (a £ - a £ )/2) / VdelM , (Al) 

which has its maximum value at the estimate ax of the sig- 
nal vector and M is the covariance matrix defined by M = 
^(a x - ax) (a £ - Sx) r }- The diagonal terms of the covariance ma- 
trix M are the squared errors Aax of the estimate ax and non- 
diagonal elements of the matrix represent correlations between er- 
rors. However, to study correlations a more convenient representa- 
tion of the correlation matrix is used - the covariance matrix nor- 
malized by its diagonal terms, My/ yfM.uM.jj ( i, j e [5, . . . , (t rec + 
l) 2 ] ). Then, the coefficient of correlation between uncertainty in a, 
and uncertainty in a, is a number between -1 and 1. 
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Figure 1. The rms of the difference between the input and output multipole maps as a function of f Iec . Filled circle, empty circle, square and triangle denote 
Cf ot for multipoles { = 2, 3, 4, 5, respectively. 
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Figure 2. Comparison of the power spectra of the WMAP map, masked by the Kp2 mask, restored by the Wiener (open circles) and PE (filled circles) filters. 
Dash dot dot and dash dot lines correspond to the averaged power spectra derived from the Wiener and PE filter reconstruction of 1() 4 simulated WMAP maps 
respectively. The smooth curve shows the WMAP team best-fitting ACDM model (£1a = 0.73 and fi m = 0.27) with running spectral index. The grey band 
indicates the cosmic variance errors (95 per cent confidence level). 
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Figure 3. Comparison of the accuracy of estimation of the aif modes using the direct inversion and PE filtering reconstruction methods, determined from 
simulations with the Kp2 mask. The abscissae give the input values of harmonic coefficients used to generate the simulated skies. The ordinates give the output 
values from the direct inversion (left column) and PE filtering (right column) methods. The coupling matrix was truncated at £ rec = 10 in the direct inversion 
method. 
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Figure 4. The same as Fig. J3J but with the coupling matrix truncated at f lcc = 30 in the direct inversion method. 
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Figure 5. Maps of the WMAP I = 2,..., 5 multipoles, computed by the PE filtering of the WMAP corrected co-added map (first row), TOH cleaned map 
(second row) and ILC map (third row) after applying the Kp2 mask. 
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Figure 6. Angular power spectrum of the WMAP data. The power spectrum from the Hinshaw et al. (2003) cut sky WMAP analysis (open circles) is compared 
with those obtained by the PE filtering of the TOH (filled diamonds), ILC (open diamonds) and corrected co-added map (filled circles), masked by the Kp2 
mask. The error bars (68% confidence level) of our WMAP power spectrum estimates are shown. They do not take into account cosmic variance. The major 
contribution to the errors comes from the correlations between the restored multipoles and between the restored and higher multipoles. The noise contribution 
is negligible. The smooth curve shows the WMAP team best-fitting ACDM model (Cl\ = 0.73 and Q. m = 0.27) with running spectral index. The grey band 
indicates the cosmic variance errors (95 per cent confidence level). 
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Figure 7. Sensitivity of the derived WMAP angular power spectrum to the input power spectrum used in the construction of the PE filter. Filled circles 
correspond to the spectrum of the multipoles obtained assuming the best fit to the WMAP data ACDM model power spectrum (indicated by the solid line). The 
open circles and filled diamonds correspond, respectively, to the same assumed power spectrum but after adding and subtracting the square root of the cosmic 
variance (indicated by the shaded band, 95 per cent confidence level). 



Figure 8. Maps of the ( = 2, . . . , 5 WMAP multipoles determined from the different frequency bands. 
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Figure 9. The angular power spectrum of the WMAP data for ST-band (filled circles), Xa-band (open circles), Q-Band (filled diamonds), V-Band (open 
diamonds) and IV-Band (triangle) maps obtained by PE filtering for the Kp2 mask. The smooth curve shows the WMAP team best-fitting ACDM model 
(Ha = 0.73 and Q m = 0.27) with running spectral index. The grey band indicates the cosmic variance errors (95 per cent confidence level). 



Figure 10. Maps of the t = 2, ... ,5 WMAP multipoles for different masks: Kp2, KpO, 20+ and 30. 
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Figure 11. The angular power spectra of the WMAP data obtained by the PE filtering of the corrected co-added Q+V+W map after masking with the Kp2 
(filled circles), KpO (open circles), 20+ (filled diamonds) and 30 (open diamonds) masks. The smooth curve shows the WMAP team best-fitting ACDM model 
(fiA = 0.73 and fi,„ = 0.27) with running spectral index. The grey band indicates the cosmic variance errors (95 per cent confidence level). 
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Figure 12. Comparison of the accuracy of estimation of the atj modes using the PE filtering reconstruction method. The abscissae give the input values of the 
harmonic coefficients used to generate the simulated skies. The ordinates give the output values determined after application of the Kp2 mask (first column), 
KpO mask (second column), 20+ mask (third column) and 30 mask (fourth column). 



Figure 13. The A^ rror maps (see text) corresponding to the sum of the errors for the restored ( = 2, . . . , 5 WMAP multipoles (centre) and for each multipole 
separately. The maps are shown in logarithmic scale. 
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Figure 14. Decomposition of the covariance matrix for the PE filtering of the WMAP data into specific contributions is shown. The total covariance ma- 
trix (upper left figure) is compared with the components induced by: correlations between restored multipoles (upper right figure), correlations of restored 
multipoles with higher order multipoles (bottom left figure) and noise correlations (bottom right figure). Only that part of the matrices corresponding to the 
range I = 2, ... ,5 are shown. The noise component is multiplied by 50 to render it visible. Maximal correlation is at most ~ 14 fiK 2 and minimal down to 
11 fiK 2 . The modes are ordered such that index i is given by i = l 2 + (. + m + 1. 



Figure 15. The multipoles map obtained by the PE filtering of the CMB maps after masking by the Kp2 mask. The CMB maps were derived from the WMAP 
data as corrected by four variants for the free-free emission. From top to bottom, the maps for following variants of the correction are shown: for Ho- template 
uncorrected for dust absorption and no free-free haze correction, for Ha template corrected for dust absorption and no free-free haze correction, for Ha 
template uncorrected for dust absorption and free-free haze correction, for Ha template corrected for dust absorption and free-free haze correction. Spectral 
indices of the thermal dust and synchrotron emission were fa = 2.2 and fa = -3.1, respectively. 



Figure 16. The t = 2 quadrupole component obtained by PE filtering of the CMB maps after masking by the Kp2 mask. The CMB maps were derived from 
the WMAP maps corrected for the free-free emission (with Ha template corrected for dust absorption and free-free haze template) for various spectral indices 
of the dust fidusi and synchrotron fa yc i r indices. 



Figure 17. Maps of the t = 2, . . . , 5 multipoles recovered using the PE filter for COBE-DMR (the first column) and WMAP data (the second column). In both 
cases the extended Galactic cut was used. 
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